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Difltusion effect on chirality selection in a two-dimensional reaction-diffusion model is studied by 
the Monte Carlo simulation. The model consists of achiral reactants A which turn into either of the 
chiral products, R or S, in a solvent of chemically inactive vacancies V. The reaction contains the 
nonlinear autocatalysis as well as recycling process, and the chiral symmetry breaking is monitored 
by an enantiomeric excess (j). 

Without dilution a strong nonlinear autocatalysis ensures chiral symmetry breaking. By dilution, 
the chiral order decreases, and the racemic state is recovered below the critical concentration Cc. 
Diffusion effectively enhances the concentration of chiral species, and Cc decreases as the diffusion 
coefRcient D increases. The relation between </> and c for a system with a finite D fits rather well 
to an interpolation formula between the diffusionless(D = 0) and homogeneous [D — oo) limits. 

PACS numbers: 05.90.-|-m, 82.40.Ck, 05.50.+q 



I. INTRODUCTION 



It has long been known that natural biomolcculcs 
choose one type of enantionicr among the two possible 
stereostructures, and the chiral symmetry is broken in 
life: Amino acids in proteins are all L-type, while sugars 
in nucleic acids are all D-type HSU. There are many 
studies on the origin of this homochirality, but the asym- 
metry caused by various proposed mechanisms turned 
out to be very minute ^ |M |y| - Frank found more than 
fifty years ago that a linear autocatalysis accompanied 
by mutual inhibition are necessary in an open system 
for the amplification of enantiomeric excess (ee) of one 
type The model with enantiomeric cross-inhibition 
is extended to polymerization models . Experimen- 
tal realization of the ee amplification is found recently in 
pyrimidyl alkanol , and explained by a nonlinear au- 
tocatalytic process in a closed system [IjI . Addition 
of a recycling process to this nonlinear system is shown 
to realize an asymptotic selection of a unique chiral state 
|l3l | . The importance of recycling process is also shown 
in nonautocatalytic but a polymerization-epimerization 
model El. 

All the above arguements assume that the reaction 
takes place homogeneously. For spatially extended 
systems, however, the inhomogcneity affects symmetry 
breaking in general. In condensed matter systems, for 
instance, elements are located densely in contact with 
each other, and interactions propagate through neigbor- 
ing elements so as to cause coherent phase transitions. 
In contrast to the equilibrium symmetry breaking in the 
condensed matter, organic molecules in the primordial 
era might have been distributed dilutely in a solution. 
If this were the case, molecular transport is necessary in 
order to bring molecules come into contact to activate 
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autocatalysis. 

In a previous letter[T5l| we briefly reported the diffusion 
effect on the chiral symmetry breaking. Here we give a 
more detailed work of the diffusion effect on the chiral- 
ity proliferation in an extended system. After introduc- 
ing the model and Monte Carlo simulation algorithm in 
§2, we consider first the chirality selection in a diffusion- 
less system, and point out some similarity of the present 
chemical reaction system to the Ising ferromagnet in §3. 
There, the effect of site dilution is also discussed. In 
§4, the diffusion effect is studied systematically and the 
chiral symmetry breaking is found to be enhanced by the 
diffusion. It is qualitatively interpreted that the diffusion 
effectively increases the concentrations of chemically ac- 
tive molecules, and consequently promotes autocatalysis. 
The idea is summarized in a simple interpolation formula 
as for the diffusion coefficient dependence of the order 
parameter and the concentration relation. The study is 
summarized in §5. 



II. MODEL AND SIMULATION ALGORITHMS 



The model we study is essentially the same with the 
one proposed in a previous letter |1!tI |. There are four 
types of molecules; an achiral reactant A, two types of 
product enantiomers R and S, and a solvent molecule in 
a diluted system, which we call a vacancy V. Molecules 
are treated as entities which occupy square lattice sites 
and can move randomly thereon. Double occupancy of a 
lattice site is forbidden. 

Various processes in the system are summarized in 
Fig. 1 . An isolated achiral reactant A turns into a chiral R 
or S molecule with a rate ko, as in Fig. 1(a). When there 
is one R(S) molecule in the nearest neighbor of a reactant 
A, the production rate of R(S) increases by an additional 
linear autocatalytic rate ki, as in Fig. 1(b). When there 
are two or more R's(S's) in the nearest neighborhood, 
the production rate of R(S) is further enhanced by an 
additional nonlinear autocatalytic rate fc2, as in Fig. 1(c). 
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FIG. 1: Various processes with transition probabilities; (a) 
non-autocatalytic, (b) linearly autocatalytic, (c) nonlinearly 
autocatalytic reactions with recycling and (d) a vacancy dif- 
fusion. X in the diffusion process (d) represents chemically 
active species, A, R or S. 



The rate of back reaction A is assumed constant, indepen- 
dent of the environment around chiral molecules, R or S 
(Fig.l(a)-(c)). By these chemical reaction processes, the 
number densities a, r and s of active molecules A, R and 
S vary, but the total number density of active molecules 
a+r+s = cis fixed constant, reflecting the assumption of 
a closed system corresponding to the Soai reaction [lol| . 
The degree of chiral symmetry breaking is represented by 
the enantiomeric excess (ee) order parameter 



(1) 



in the list has the same transition probablity Pi to per- 
form the reaction or diffusion process i, shown in Fig.l. 
Assume that there are Ui lattice sites in the correspond- 
ing list of the process i. Then, the total probability of 
events in a unit time is obtained as Ptot = X^i ^i'^i- 
other words, in a time interval dt = 1/Ptot, one of the 
process takes place on average. By generating a quasi- 
random number x which is uniformly distributed between 
and Ptot 7 the process j which satisfies the condition 



< a; < ^ P,5 



(2) 



is to be performed. The lattice site which performs the 
process j is listed at the fc-th position in the list of Uj 
sites, where y = x — Pi'>^i ^J^id k is the integer part 

of y/Pj -t- 1 as 



k = 



+ 1 



(3) 



After performing the process j on the site in the A:-th 
position in the list, the time is advanced by dt and the 
lists are updated. We simulate systems with sizes 50^ 
and 100^, but in the following we show results obtained 
for the larger system with a size 100^. 



III. WITHOUT DIFFUSION 

From our previous works on a homogeneous system 
ITsI the nonlinear autocatalysis with a finte ^2 is found 
essential for the chiral symmetry breaking, whereas the 
linear autocatalysis ki is inefficient. Therefore, fci is set 
absent hereafter: fci = 0. In order to find out the effect 
of diffusion, we consider in this section the developement 
of chirality in space without the material transport. The 
effect of diffusion is studied in the next section. Here, 
the order of chiral symmetry breaking is studied as a 
function of the strength of the nonlinear autocatalysis or 
of the concentration of active molecular species. 



A vacancy V is inactive as for the chemical reaction, 
but plays a vital role in the diffusion process. When 
chemically active species. A, R and S, are neighboring 
to V, they exchange their sites with V randomly with 
a life time td on average, as in Fig. 1(d). The diffusion 
coefficient is then D ~ /{Atd) with a being the lat- 
tice constant and set to be unity hereafter. Previously 
we simulated oirly the case with D = 1/4, because the 
Metropolis algorithm was used in the Monte Carlo sim- 
ulation [l5l |. In the present paper we want a system- 
atic study of diffusion effect on the chirality selection, 
and therefore an arbitrary variation of the diffusion co- 
efficient D is necessary. For that purpose, we adopt a 
waiting time algorithm for the Monte Carlo simulation. 

The simulation runs as follows: All the positions which 
can undertake a process i are tabulated in a list; every site 



A. Chiral Symmetry Breaking as a function of 
Nonlinear Autocatalysis 

First, we consider the A:2-dependence of the chiral sym- 
metry breaking, when the whole lattice sites are occupied 
by chemically active species; c = 1. When the nonlinear 
autocatalytic effect is strong enough such as k2 — lOOfep 
with A = fco, achiral reactant A quickly reduces to chi- 
ral species R or S. Therefore, even by starting from an 
achiral initial configuration with only achiral reactants 
A, domains of R and S are formed in the early stage, 
as shown at a time k[)t = 50 in Fig^ Then, the do- 
main competition sets in, and those domains surrounded 
by another type shrink as shown in the middle column in 
Fig. 2. Eventually, competition ends up by the dominance 
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FIG. 2: Configuration evolution of (a) R and (b) S molecules 
in a difTusionless system. Tmes are kot = 50, 250, 1750 from 
left to right. Parameters are ko = 
c = 1, and the system size is 100^ 
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FIG. 3: Enantiomeric excess |(/)| versus ki/ko at the full oc- 
cupation c = 1 in a semi-logarithmic plot. A curve represents 
the temperature dependence of the magnetization 1201 of a 
square Ising model, where a temperature is related to the 
chemical reaction rates as fe2/feo = exp(4J/fcBT) — 1. 



of one enantiomeric type, S in the case shown in the right 
column of Fig. 2. In some cases the system decouples into 
two inter-penetrating domains of R and S, both of which 
in effect extend to infinity due to the periodic boundary 
conditions in the x and y directions: Then, it takes an 
enormous time until the minor domain extinguishes and 
the system accomplishes the final equilibrium state. 

In order to obtain the equilibrium value of the order 
parameter (/> in a reasonable simulation time, therefore, 
we start simulations from an initial configuration such 
that the whole lattice sites are occupied by R species. 
Equilibrium states are readily realized from this initial 



FIG. 4: Analogies of state transition in the reaction model 
(above) and the Ising ferromagnet (below), (a) n = 2, and 
(b) n = 3. 



configuration, and the equiblirium value of the order pa- 
rameter (j) thus obtained is shown in Fig|31 at various 
strength of nonlinear autocatalysis fc2/fco- Since the non- 
autocatalytic reaction fcp drives the system to the racemic 
state, the value of the order parameter <f> decreases as fc2 
decreases, and vanishes at a critical value of k2/kQ around 
10: The system recovers the chiral symmetry, getting into 
the racemic state. 

The qualitative features of domain competition shown 
in Fig. 2 and the emergence of symmetry breaking at large 
fc2 as shown in Fig. 3 suggest us some resemblance with 
the phase transition in an Ising ferromagnet. In fact, 
one may find some analogy of our reaction system to 
an Ising ferromagnet: In the latter, each lattice site is 
occupied by an Ising spin, which can be in one of the two 
states; up or down. When a nearest neighboring spin pair 
have the same orientation, the energy is low as given by 
— J < 0, whereare when a pair is in opposite orientations, 
the energy is high as J > 0: The energy difference is 
set to be 2J > 0. At a temperature T, a spin flips its 
orientation with a transition probability which depends 
on the associated energy change. Let the spin to flip have 
n nearest neighboring spins with the same orientation. 
Then, it will have 4 — n spins of the same orientation 
after the flip, and we denote the transition probability as 
W{n A ^ n). In order to ensure thermal equilibrium 
asymptotically, it is known to be sufficient if the detailed 
balance condition between the transition probabilities of 
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one process and its reverse 
W{n 4- n) 



W{4 



n) 



cxp(-A£:/fcBT) 



(4) 



is satisfied, where AE = 4(n — 2) J represents an energy 
change. For example, when a site is surrounded by two 
up spins and two down spins, n ~ 2 and the transition 
probabiUty of an up spin to flip downward is equal to its 
reverse process, as shown in the lower column of Fig. 4(a). 
Then, in equilibrium the probabilty P+ (2, 2) of an up spin 
surrounded by two up spins and two down ones is equal 
to that of a down spin in the same enviornment P_ (2, 2): 
P+(2, 2) = P_(2, 2). We associate, for instance, R and S 
molecules to the up and down Ising spins, respectively. 
Then, the transitions for n = 2 in the Ising model corre- 
spond to transitions between an R and S molecule sur- 
rounded in the neighborhood by two R's and two S's ( 
upper column of Fig. 4(a)). Since the alternation of R and 
S is mediated by an achiral reactant A, the equilibrium 
probabilities of R and S in the configuration Fig. 4(a) arc 
related as APi? (2, 2) = (fco + fc2)PA(2, 2) = APs(2,2), and 
thus Pr(2, 2) = ^5(2, 2), consistent with our assignment 
of R and S molecules to up and down spins. When a 
site is surrounded by three R's and one S molecule as in 
the upper column of Fig. 4(b), the equilibrium probabil- 
ity that the site is occupied by an R molecule Pr{3, 1) is 
higher than the one by an S molecule ^^(S, 1). They are 
related as APr(3, 1) = {ko + k2)PA{3, 1) and koPA{3, 1) = 
APs(3, 1), and thus Pr{3, !) = (! + fc2/fco)Ps(3, 1): At a 
lattice site surrounded by three R's and an S, the proba- 
bility of finding an R molecule is higher than that of find- 
ing an S molecules by a factor (1 + k2/ko). For the Ising 
system, the corresponding situation is the spin on the site 
surrounded by three up's and one down; n = 3. The re- 
lation between probabilities that a site is occupied by an 
up spin and by a down is e-4"'/'=^«^P+(3, 1) = P_(3, 1). 
Thus, we may assume that the ratio (1 + fc2/fco) cor- 
responds in the Ising model to the ratio of the transi- 
tion probabilities exp(4J//cBT'). (See the lower column 
of Fig.4(b)). In other word, the nonlinear autocatalytic 
rate ^2 in our chiral model might correspond to the cou- 
pling constant of the two-dimensional Ising model as 



k2/ko 



1. 



(5) 



Since the numbers of R and S molecules correspond to 
those of up and down spins in the Ising fcrromagnct, the 
chiral order parameter (f> corresponds to the spontaneous 
magnetization M of the Ising model. The latter is exactly 
known to behave as [13 M = [1 - sinh(2 J/A:i3r)-4]i/8. 
By using the relation Q , the magnetization M is plotted 
in Fig. 3 by a continuous curve. 

The transition temperature Tc of the Ising model on 
the square lattice is exactly known to be at fc^Pc/SJ — 
1/ ln(l-|--y2) ~ 1.135. Thus in the present reaction model 
the chiral symmetry breaking is expected to take place at 
^2 > ^2c = 2(-y2-|-l)fco ~ 4.83fco. However, the order pa- 
rameter (j) obtained by simulations decreases faster than 



the Ising magnetization as k2 decreases: For the chiral 
symmetry breaking a larger value of the nonlinear auto- 
catalysis k2 ~ lO/co is required, as shown in FigEl This 
is because of the approximate character of the relation 
between the Ising coupling constant J and the nonlin- 
ear autocatalytic rate k2. The production rate in our 
model is independent of the number of neighboring R(S) 
species as far as it is more than 2, whereas transition 
rate to the up spin state in an Ising model increases 
as the neighboring number of up spins increases: This 
means that our coupling constant k2 is effectively weaker 
than the Ising coupling parameter. Also the existence 
of an achiral species is not included in the Ising model. 
The decomposition process from a chiral species R or S 
back to an achiral reactant A induces fluctuation which 
breaks ordering, and a chirally ordered state in our re- 
action model is less stable than the ferromagnetic state 
in the Ising ferromagnet: The strong nonlinearity is nec- 
essary for the symmetry breaking. In fact, the critical 
value for the chiral symmetry breaking k2c depends on 
the decomposition rate A, though the ^2 versus J relation 
in eq.(5) is independent of A. From simulations with a 
large A = IO/cq- the order is found to be easily broken at 
k2c ~ 23fcQ. whereas for small A = O.I/cq the order holds 
until k2c ~ 9fco- 



B. Dilution Effect on Chiral Symmetry Breaking 

We now keep the nonlinear autocatalitic effect as 
strong as k2/ko = 100, and consider the dilution effect 
(c < 1) on the chiral symmetry breaking. If there is no 
material transport, the behavior of the system depends 
strongly on the arrangement of vacant sites. There- 
fore, in order to obtain the order parameter (j) we av- 
erage 9 samples at concentrations lower than c = 0.7, 
whereas 3 samples seem sufficient at higher concentra- 
tions. As the system is diluted, the probability that ac- 
tive chemical species are contiguous to each other de- 
creases. In the present reaction model, however, a reac- 
tant A should be surrounded by chiral products R or S to 
sustain the autocatalysis. Therefore, the chiral symmetry 
breaking is possible only at concentrations at least larger 
than the site percolation threshold of the square lattice 
Cp = 0.592746 Below Cp, the information of one 

site is unable to propagate through the whole system, and 
thus the system cannot be in a single state. The order pa- 
rameter (j) actually decreases as the system is diluted, as 
shown in Fig. 5a, and the system becomes racemic at con- 
centrations lower than the critical value Cc{D = 0). By 
assuming a critical concentration at Cc{D = 0) « 0.639, 
the order parameter (j) shows a power-law dependence 
4> o: \/ c — Cc{D = 0), as shown in Fig. 5(b). The critical 
concentration Cc{D — 0) without diffusion is close to but 
larger than Cp. 

A rather large discrepancy between Cc(P = 0) and 
Cp might be understood by using the similarity of the 
present diluted reaction model to a site-diluted Ising for- 
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romagnct, where lattice sites are diluted substitution- 
ally by nonmagnetic impurities ^ site-diluted 
Ising model, the transition temperature is known to de- 
crease as the concentration of nonmagnetic impurity in- 
creases, and the phase transition vanishes at the criti- 
cal concentration close to the percolation threshold Cp 
jl9l |: The transition temperature is zero at the critical 
concentration. If the relation (5) holds, however, our re- 
action model with ^2 = lOOfco corresponds to an Ising 
ferromagnet at a finite temperature ksT /2J = 0.433 or 
T/Tj, = 0.38. A thermal fluctuation favors the recoveray 
of the symmetry even if a percolating cluster exists, and 
this may be the reason that the critical concentration 
Cc{D = 0) is larger than the percolation threshold Cp. 
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FIG. 5: Enantiomeric excess \(f>\ versus c in a difTusuioness 
system with a strong nonlinear autocatalysis ^2/^0 = 100. 



IV. DIFFUSION ENHANCEMENT OF THE 
CHIRALITY SELECTION 

We now study the diffusion effect on the chirality se- 
lection in a diluted system. The transport of material 
in general increases the chance for an achiral reactant A 
to meet chiral products R and S in their lifetimes, and 
may enhance the autocatalytic effect. When the diffu- 
sion is extremely fast, the system becomes homogeneous: 
This case will be studied in the following subsection an- 
alytically. With a finite diffusion coefficient, numerical 
simulation is the only possible method to study, and the 
results so obtained will be analyzed by the interpolation 
formula in the second subsection. 



A. Homogeneous Case with D = 00 

When the diffusion is extremely fast, the system be- 
comes homogeneous. Then, the simple mean field ap- 
proximation leads to the rate equations for the homoge- 
neous concentrations; 

r — k(r)a — Xr 

s — k{s)a — As, (6) 

with the rate constant k{r) = ko + 4fcir -I- 6fc2r^. Here 
in this subsection we restore the rate fci of the linear au- 
tocatalytic reaction, in order to be general. The factors 
4r and 6r^ represent the probabilities that at least one 
or two of the four neighboring sites are occupied by R 
species, respectively, in the sense of mean field approx- 
imation. Since the system is assumed to be closed, the 
total concentration of chemically active species satisfies 
the concervation law; a + r + s = c =const. Then, the 
total concentration of chiral products u = r + s is shown 
to varies as 

u = [2ko + Akiu + 3fc2U^(l + (^^)](c - u) - Aw, (7) 
whereas the order parameter evolves as 

0= ^^</)[-2fco + 3A:2u'(l-0')]. (8) 
u 

Since < u < c, eq.(8) is similar to the time-dependent 
Landau equation with nontrivial fixed points if > 
2ko/{3k2). The other point to remark is that eq.(8) does 
not depend explicitly on the linear autocatalytic rate fci ; 
(j) depends on ki only implicitly through u. Thus the 
linear autocatalysis alone cannot break the chiral sym- 
metry. 

The equilibrium values of u and are calculated from 
the algebraic equations given by (/) = u = 0. The sponta- 
neous chiral symmetry breaking takes place at a critical 
concentration Cc{D = 00) when the diffusion is infinitely 
fast: The order parameter (f> of a, homogeneous system 
takes a finite value for Cc{D = 00) < c < 1. For fci = 
the order parameter (f) is explicitly written as a function 
of the concentration c as 

and the critical concentration is calculated to be Cc {D = 
00) = (4 A/fco)/^fco/24fc2. It takes the value Cc{D = 
00) = 0.1020 for the parameter values ki = 0, k2 — 
lOOfco, A = fco. 

B. Finite Diffusion Constant 

The emergence of chiral order at a finite vacancy diffu- 
sion is studied by the Monte Carlo simulation. It is found 
that the symmetry breaking takes place rather rapidly 
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FIG. 6: (a) Enantiomeric excess \<j)\ versus c at various values 
of the diffusion coefficient D/ko. The reaction parameters are 
k2/ko = 100, fci = and A = fco- A continuous curve rep- 
resents the relation (9) for a homogeneous system. Dashed 
curves represent interpolation (10) at various diffusion coef- 
ficients, (b) Diffusion coefficient dependence of the concen- 
tration c where the order parameter takes the value (j> = 0.6. 
A curve represents the interpolation (10) with a parameter 
A = 0.235. The dotted line represents the homogeneous case 
with c(0.6; D ^ oo) = 0.118. 



even if we start from a completely achiral initial con- 
figuration with all the chemically active molecules be- 
ing achiral reactant A as a = c: The order parameter 
approaches quickly the asymptotic value, though which 
enantiomer dominates over the other is determined ran- 
domly in each simulation run. Therefore, we initiated 
the simulation from the complete achiral state with only 
A and V species. 

The steady-state value of ee or \(f)\ at various values of 
diffusion coefficient D is shown as a function of the con- 
centration of active species c by symbols in Fig|B{a). A 
continuous curve represents the relation (9) for a homo- 
geneous system with D = oo. The most striking feature 
is that (j) — c relations at different D's look similar but 
with shifts to the lower concentration as D increases. 
The D-dependence of the shift might be considered in 
the following way. The chiral products R or S once pro- 
duced diffuse around in a diffusion area of about D/X 
sites during the life time A. Therefore, in a system where 
diffusion is allowed the effective density of chiral cata- 
lysts is enhanced by a factor 1 -I- AD/X in comparison to 
the system without diffusion. A proportionality constant 



A will be determined later by fitting to the simulation 
results. With a very large diffusion D, however, there 
will be overlappings of various diffusion areas, and the 
diffusive enhancement may saturate. Thus, even with an 
infinitely fast diffusion, a finite concentration of active 
species is necessary, as is given by eq.(9). 

The diffusive enhancement of the effective concentra- 
tion means that in order to reach a fixed value of the 
order parameter 0, the required concentration of active 
species c(^; D) for a system with a finite diffusion coef- 
ficient D can be less than c((/); D = 0) for a diffusionless 
system. Here we propose a simple formula for c(^; D) 
interpolating between the diffusionless limit c{(f>; D = 0) 
and the homogeneous limit c{(f>; D ~ oo) as 



b;D) 



c(0;O)+c(^;oo)(Ai^/A) 
1 + [AD/X) 



(10) 



where c((/); 0) is obtained by the simulation (Fig. 5(a)), 
whereas c{(j)] oo) is obtained by solving eq.(9) as 

4fco 



To determine the constant A, we fit the formula (10) with 
simulation data at = 0.6. The reason for the choice of a 
non-zero value of is that the value of c{(j}; D) is expected 
to be more accurate away from the transition point, 4> = 
0. The fit gives the value A « 0.235, and agrees well at 
small D but deviates systematically at large D, as shown 
in FigEtb). 

With this value of A, we interpolate the whole (j) — c 
curve for arbitrary D by eq.(lO), as is drawn by broken 
curves in Fig. 6 (a). The trend of the down shifting in 
the concentration as the diffusion increases is well repro- 
duced, even though the fitting is not perfect, especially 
at a large diffusion coefficient D. There, we have so far 
no plausible theory on the Z3-dependence of c{(j)\ D). 



V. SUMMARY AND DISCUSSIONS 

We study diffusion effect on the chirality selection dur- 
ing a chemical production of chiral molecules, R or S, 
from an achiral substrate A in a solution by means of a 
square lattice model. In addition to the usual reactions 
of the production and the decompsition of chiral species, 
autocatalytic processes are incorporated such that an 
achiral substrate surrounded by many enantiomers of one 
type increases the turn-over rate to that type. Specifi- 
cally, the nonlinear autocatalysis which takes place when 
the substrate is surrounded by two or more enantiomers 
of one type is found to be essential in the chiral symmetry 
breaking. 

When all the lattice sites are ocuupied by chemically 
active species, A, R or S, the spontaneous symmetry 
breaking has some similarity to the phase transition of 
an Ising ferromagnet. The strength of the nonlinear au- 
tocatalysis k2 in our chemical model is approximately re- 
lated to the Ising coupling constant. The chemical model 
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at very large ^2 evolves through the domain competition, 
similar to the dynamics of the Ising model at low temper- 
atures. As k2 decreases, the chiral order extinguishes, but 
earlier than the expectation from the Ising model. The 
discrepancy may be due to the larger degrees of freedom; 
three in the chemical model whereas two in the Ising 
model. 

Inclusion of vacancies or dilution suppresses the degree 
of the chiral symmetry breaking or the enantiomeric ex- 
cess (f), and below the critical concentration of chemically 
active species, the system is in a racemic state with a 
chiral symmetry. In a diffusionless system {D = 0) with 
a strong nonautocatalysis, the system is similar to the 
site-diluted Ising model, and the critical concentration is 
close to but larger than the percolation threshold Cp. 



Diffusion enhances the chiral order 0, since the migra- 
tion of chiral products enhances the chance of molecular 
collision and the autocatalytic process as well. In the 
limit of an infinitely fast diffusion [D = oo), the system 
behaves as if it is homogeneous and the rate equation in 
the mean field approximation describes the chiral sym- 
metry breaking: The steady state solution is obtained 
analytically in this limit. At a finite diffusion coefficient 
D, the concentration dependence of (f) is shifted from the 
diffusionless limit {D = 0) to the homogeneous limit 
{D = oo). The behavior is interpreted as the diffusive 
enhancement of concentrations of chiral species, and is 
summarized in an interpolation formula of c-dependence 
of (j) or c(0; D) between the two limits. 
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